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Abstract. We calculate the Landau interaction function /(k, k') for the two-dimensional t-t' Hubbard 
model on the square lattice using second and higher order perturbation theory. Within the Landau-Fermi 
liquid framework we discuss the behavior of spin and charge susceptibilities as function of the onsite 
interaction and band filling. In particular we analyze the role of elastic umklapp processes as driving force 
for the anisotropic reduction of the compressibility on parts of the Fermi surface. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 74.72.Jt Other cuprates 



1 Introduction 

The Hubbard model remains one of the most studied the- 
oretical models in the efforts to understand the high-tem- 
perature superconducting cuprates. Despite their excep- 
tionally high transition temperatures into the supercon- 
ducting state, much of the interest in these materials arises 
from their anomalous normal state properties over a large 
part of the temperature-doping phase diagram. These ef- 
fects occur unequivocally and most strongly in the under- 
doped samples with electron densities close to half band 
filling where the system is a stoichiometric antiferromag- 
nctic Mott insulator. In the undcrdopcd regime, the tran- 
sition into the superconducting state occurs out of the 
so-called pseudogap phase, which shows significant devia- 
tions from conventional metallic Fermi-liquid behavior Jj], 
||. This is in contrast with sufficiently overdoped sam- 
ples where many observations point toward a conventional 
transition from a Fermi-liquid into the superconducting 
state There have been experimental and theoretical 
studies suggesting that a quantum critical point hidden 
underneath the superconducting dome j|,|j||6| or the onset 
of Fermi surface truncation fj|[§J| give rise to these phe- 
nomena. On the other hand one may try to take a more 
conservative route and attempt to understand the experi- 
mental findings within a conventional Fermi-liquid frame- 
work. Typically, rather specific assumptions on the Lan- 
dau interaction function have to be madcjio) in order to 
obtain agreement with the experimental data. Thus it ap- 
pears useful to simply calculate within perturbation the- 
ory the Fermi-liquid properties of the 2D Hubbard model 
and to see whether these assumptions can be justified 
within this transparent theoretical framework. 
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Another motivation for this work is the question whether 
as the density is increased towards half filling some kind 
of drastic change of the low energy properties, e.g. corre- 
sponding to Fermi surface truncation as indicated by one- 
loop renormalization group calculations is foreshad- 
owed already in a perturbative Fermi liquid picture. Such 
a partial truncation of the Fermi surface has been argued 
to be seen in angular resolved photoemission experiments 
0. It is a promising candidate to understand the many 
anomalies of the underdoped normal state of the high- 
T c cuprates, as it naturally combines insulator-like fea- 
tures such as the pseudogap, the observed c-axis resistivity 
and the strongly reduced supcrfluid weight, with residual 
metallic properties as evidenced by the in-plane transport 
and superconductivity itself. Since both the renormaliza- 
tion group approaches and the present Fermi liquid analy- 
sis are perturbative techniques, we will not obtain an accu- 
rate description of possible phases with a truncated Fermi 
surface, but we will carefully look for insulator-like tenden- 
cies in these calculations as the interaction gets stronger. 
In particular we analyze the influence of umklapp pro- 
cesses between k-space regions in the vicinity of the sad- 
dle points at (ir,0) and (0, 7f). Elastic umklapp scattering 
with momentum transfer of (jr, 7r) across the Fermi surface 
is allowed when the Fermi surface extends to the Brillouin 
zone boundary for a certain density range close to half 
filling. 

A Fermi liquid analysis of the 2D Hubbard model close 
to half filling has been presented by Fuseya et al.pi|. 
These authors studied the case with only nearest neigh- 
bor hopping, t, and found a suppression of the uniform 
spin susceptibility for sufficiently strong interaction close 
to half-filling, which maybe interpreted as indicating the 
opening of a spin gap. Unfortunately this promising indi- 
cation turns out to be an artifact of the calculation scheme 



used by Fuseya ct al. 11 , as will be discussed below. From 



2 



Frigeri, Honerkamp, and Rice: Landau-Fermi liquid analysis of the 2D t-t' Hubbard model 



renormalization group studies @| we know that nonzero 
values of the next nearest neighbor hopping, t' , can change 
the weak coupling physics significantly. Therefore, in this 
paper we extend the analysis of the Landau interaction 
function and related quantities to the case of t' ^ 0. Par- 
ticular attention is paid to the anisotropy of the interac- 
tion function for different parts of the Fermi surface. 

In the following we first apply a second order pertur- 
bative method and discuss the k-space integrated and k- 
space resolved Fermi liquid properties of the 2D t-t' Hub- 
bard model. Next we describe how ladder summations in 
the particle-particle and particle-hole channels, which are 
necessary to remove unphysical divergences, modify the 
tendencies observed in the second order approach. Finally 
we conclude and compare our results to other approaches 
and experimental observations for the high-T c cuprates. 



which expresses the fact that r has to be antisymmet- 
ric with respect to the exchange of two particles with the 
same spin orientation. Then the Landau interaction func- 
tion can be obtained via (|) by taking the limit (|) in the 
perturbation expansion for V SS /(K, K',K + Q). The first 
and second order graphs contributing are shown in Fig. 
[j]. First the Landau function will be evaluated to second 
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2 Second order calculation of the 
quasiparticle interaction 

The kinetic energy for the t-t' Hubbard model on the 2D 
square lattice is given by 
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ek = — 2i(cos k x + cos k y ) + At' cos k x cos k y 



(I) 



where t and t' denote the amplitudes for hopping be- 
tween nearest and next-nearest neighbors, respectively. 
The short range Coulomb interaction between the elec- 
trons is taken into account by the usual onsite repulsion 



H, 



u 22 n > " 



(2) 



The particle number N e is controlled by a chemical po- 
tential (i. 

The Landau interaction function / s . s /(kp, k' F ) can be 
obtained from the two-particle vertex r ss > (K, K', K + 
Q). Here K = (k, v) is a short notation for wavevec- 
tors and frequencies. F is a function of the interaction 
strength, U. To first order in [/, r is wavevector indepen- 
dent, but higher orders lead to a pronounced wavevector 
dependence. K, K' and spin indices s,s' characterize the 
two initial particles and K + Q, K' — Q and s,s' their final 
states, and Q = (q, uj) denotes wavevector and frequency 
transfer in the scattering process. The Landau quasiparti- 
cle interaction /(k F ,k F ) corresponds to the vertex func- 
tion in the limit Q — > with £ — > with K = (k F , 0) and 
K' = (k' F ,0). More precisely" 

/ s , s ,(k F ,k' F ) = z k z k , lim limr s , s ,(K,K',K + Q), (3) 

UJ— >0 q^O 

where k F , k' F are two wavevectors on the Fermi surface 
and Zfc is the residue of the pole of the interacting one- 
particle Green's function G(k F , v), i.e. z k = [l—SS/Sv)^ 1 
with self-energy £(lt, v). The two-particle vertex is de- 
termined by the coupling function V^ S /(K,K',K + Q) = 
y(K,s;K',s';K + Q,s) via 
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Fig. 1. First and second order graphs contributing to F. We 
use the notations j) = { , } and f = {', }. 



order in U. Corrections from the z k factors are at least 
second order in U, therefore setting z k = 1 is fully consis- 
tent up to second order. The same argument justifies the 
neglect of effective mass corrections. Finally the second 
order quasiparticle interaction function reads 



/ TT (k F ,k F ) = / U (k F ,k F ) 



(5) 



/ Ti (k F ,k F ) - fu(k F ,k' F ) 

= U + U 2 (xpp(^f + k' F ) + xph (k F - k' F )), 

(6) 
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Here n p = 1/(1 + cxp(e p — fi) /T) is the Fermi distribution 
at temperature T . 

For the analysis of the physical properties of the in- 
teracting Fermi liquid it is convenient to introduce sym- 
metric (f s ) and antisymmetric parts (f a ) of the Landau 
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quasiparticle interaction by 

1 



f' a - o(/TT ± /it)- 



(9) 



/ s (ki?,k^) and f a (kp, k' F ) are computed for a fixed 
value of the chemical potential fi. kp and are parame- 
trized by the two angles 6 and 0' ', measured from the sym- 
metry axis FX, as shown in the inset of Fig. g. Fig.|] 
shows the results obtained for the parameters t' Jt = 0.3, 
n/t' = -3.83, T/t = 0.002 and U/t = 1. Note that for 
this choice of the chemical potential fi/t' > —4, the Fermi 
surface consists of four arcs which lead to the breaks in 
the curves. The most prominent feature is the divergence 
along a diagonal line with it < 9' < 3tt/2 in both func- 
tions. This reflects simply the cooper divergence in second 
order perturbation theory and it is removed by summing 
the ladder graphs in the particle-particle channel as dis- 
cussed below. 
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3 Total charge and spin susceptibility for the 
almost isotropic case 

In a Fermi liquid the static uniform charge and spin sus- 
ceptibilities X c and X s involve only excitations in the 
neighborhood of the Fermi surface. Therefore they can 
be calculated as Fermi surface integrals over k-space local 
quantities Xc{^-f) and XsO^-f) which define local suscepti- 
bilities for Fermi surface points described by the wavevec- 
tor kp, 







Fig. 2. Numerical results for the symmetric f s (9,0') and an- 
tisymmetric part f a (8,9') of the Landau interaction obtained 
by the second order perturbation approximation for the pa- 
rameters t' = 0.3t, fi = -3.83t', U = It and T = 0.002t. To 
the second order / a,3 (kF,— kf) show unphysical divergences 
(sharp diagonal features for ir < 6' < 37r/2). 



P 2 X, 



1 dN P 



X s 1 dM 
fi 2 B p B Q dH 



dl F XcQlf), 
dlF Xs{^f) 



(10) 
(11) 



Here, J dip denotes the line integral along the Fermi sur- 
face and Q is the total volume. N e denotes the number 
of electrons and p — N e /f2 denotes their density. [Ib is 
the Bohr magneton. The k-space local susceptibilities are 
solutions of two inhomogeneous linear equations, 



L s (xc) = 2, L a ( Xs ) = 2, 



(12) 



which express via the two linear integral operators L s 



L s ' a ( X ) 
(2tt)2 
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(13) 



how the local susceptibilities are renormalized by the Lan- 
dau /-function |L3]]. Here, vp = Vk F denotes the magni- 
tude of the Fermi velocity Vk F = 9ek/9k|k=k F of a quasi- 
particle with wavevector k. For an isotropic Fermi surface 
and isotropic interactions, Xc(k) and Xs(k) are obviously 
two constants independent of k. This allows one to ex- 
press X c and X s using the conventional Landau parame- 
ters F^ a , which are defined as zeroth order coefficients in 
the expansion of the Landau interaction functions f s ' a in 
normalized Legendre polynomials. In the anisotropic case 
the Landau interaction functions can be expanded in a 
similar way, 



r a (kp,k' F ) = -L E K:Mkp)^(k> F ) , (u) 

1 V 



71=0 7 



where ijjo(kp) = 1, tp n are orthogonal functions on the 
Fermi surface and Np is the density of states at the Fermi 
surface. Note that in general the off-diagonal terms F*'*, — 
F^} a n do not vanish when both ip n and tp' n belong to the 
same irreducible representation of the crystal symmetry 
group. 

Initially, let us assume that the isotropic components 

of ,r> a , 
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(15) 
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are dominant. This assumption was also made in Ref. 
In this case the compressibility ([n]) simplifies to 
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P 2 X C 



N F -F° J dl' F xM 
N F - F°p 2 X c . 



(16) 



Thus we arrive at the well-known result for an isotropic 
Fermi liquid, 



p 2 x c (n) = 



N F 



(17) 



Using the same simplification we obtain for the spin sus- 
ceptibility 



Mb 



N f 



l + F Q (/i)- 



(18) 



Fuseya et al. evaluated charge and spin susceptibili- 
ties for the Hubbard model with t' = assuming that the 
isotropic components of f B ' a were dominant. Their results 
showed a reduction of the uniform spin susceptibility for 
U > 2t at density near to half-filling. It turns out that this 
effect is tied the van Hove band filling where the density 
of states at the Fermi surface diverges due to van Hove 
singularities at the saddle points (it, 0) and (0, ir). In the 
special case t' = the reduction of the spin susceptibil- 
ity does not have a simple physical interpretation. In this 
case the van Hove filling coincides with half-filling, where 
the system is believed to be an antiferromagnetic insula- 
tor and the Fermi liquid approach must break down at low 
energy scales. But for finite next-nearest neighbor hopping 
t' the van Hove filling does not coincide anymore with 
the half-filled or perfectly nested case and in the vicinity 
of this filling the system may be in a paramagnetic state. 
In this case a reduction of the uniform spin susceptibil- 
ity for U > 2t could be interpreted as an indication for 
an opening of a spin gap in the Hubbard model close to 
half filling. However, our results obtained by solving the 
full integral equations (|l2|), without the assumption of a 
dominance of the n = terms, unfortunately remove this 
promising finding. This is shown in Fig. where the data 
assuming that the isotropic terms are dominant (left pan- 
els) are compared with those obtained by solving ([nf |l2|) 
(right panels) . In particular in the two figures to the right 
the spin susceptibility does not decrease in any significant 
way as the van Hove filling (/i = for t' = and p/t' = —4 
for t' — 0.25t) is approached. For the charge susceptibility 
the qualitative behavior of the results obtained by solving 
( |Tc| , p^l) is compatible with the one obtained by Fuseya, 
even if a quantitative difference is present, see Fig. 

Thus the conclusion of Fuseya et al.jyj that indica- 
tions for a spin gap in the 2D Hubbard model can already 
be observed in second order perturbation theory seems 
to be an artifact of their approximation. The solution of 
the full integral equation reveals that the assumption of a 
dominant I = Landau parameter is not always justified. 
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Fig. 3. Comparison between the uniform spin susceptibility 
X a as function of the chemical potential fi obtained from the 
approximate relations (^|) (plots on the left) and from the so- 
lution of the full integral equations ( ^j|Jlll ) (plots on the right), 
for different values of interaction U: * U/t = 1, o U/t = 2, \7 
U/t = 3, A U/t = 4. For all plots T/t = 0.002. The suscep- 
tibility for the non-interacting system (U/t = 0) is shown as 
dashed line. 
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Fig. 4. Same as in Fig. ^, but for the charge susceptibility X c 



4 Local spin and charge susceptibilities in 
k-space 

Since the Landau interaction functions are very anisotropic, 
it is interesting to analyze how the local values of the static 
susceptibilities change with the interaction U. As men- 
tioned in the introduction, we are particularly interested 
in the question whether there are tendencies towards a 
Fermi surface truncation, i.e. the formation of incompress- 
ible k-space regions at the saddle points, as suggested by 
one-loop renormalization group studies, already within a 
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low order Fermi liquid treatment. For the non-interacting 
system, XcQ^-f) and x s (kp) are proportional to the inverse 
velocity 1 / Vk F ■ Since Wk F vanishes at the saddle points, a 
divergence appears in the local spin and charge suscep- 
tibilities around these points at the van Hove filling - a 
behavior opposite to incompressibility. The results for the 
interacting system obtained with the second order Lan- 
dau function show a general reduction of the local values 
of the susceptibilities with increasing interaction strength 
U. Fig. |]shows XsO<-f) and XcO<-f) when is close to the 



points on the Fermi surface (e.g. 9 « ir/8 and 9 
while for Xc(Kf) the situation is reversed. 
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Fig. 5. k-space local values of Xs(kf) and Xc(^-f) for \s.f near a 
saddle point for different values of U : U/t = 1 thick line, U/t = 
2 thick dashed line, U/t = 3 thin line, U/t — 4 thin dashed 
line. The inset shows how the angle 9 is used to parameterize 
the points of the Fermi surface (kp(#)). Results computed for 
t' = 0.3, n/t' = -3.99, and T/t = 0.002. 



saddle point. The inset in the spin susceptibility graph in- 
dicates how the angle 9 is used to parameterize the points 
of the Fermi surface, kp(0). Particularly interesting is the 
fact that near the saddle point the reduction of the local 
spin susceptibility is pronounced. The charge susceptibil- 
ity shows a similar behavior, but in a less marked way. In 
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Fig. 6. Relative variations of the k-space local spin and charge 
susceptibilities (Xs(U) and Xc(U)) near the van Hove filling for 
three different values of the angle 9 as function of the interac- 
tion U. 9 = Q m in thick line, 6 w 7r/8 thin dashed line, 6 — 7r/4 
thin line. T/t=0.002. 



fact the comparison between the two graphs of the Fig. 
H shows that the relative value of Xs(^-f) decreases more 
rapidly near the saddle points (9 = 9 m i n ) than at other 



Fig. 7. Left plots: relative variation of the local spin and charge 
susceptibilities (xs{U), Xa(U)) for three different values of the 
angle 9 as function of the interaction U. 9 = 9 m i n thick line, 
9 « tt/8 thin dashed line, 9 = vr/4 thin line. T/t = 0.002. 
The band filling is larger than the van Hove filling and elastic 
umklapp processes at the Fermi surface are possible. Right 
plots: results obtained excluding all umklapp processes from 
the Landau interaction functions. 



The behavior of the susceptibilities for values of [i fur- 
ther away from the van Hove filling is noteworthy, too. 
In particular the qualitative difference of the susceptibil- 
ities between the case /i > —At 1 (density larger than van 
Hove filling) and /i < —At 1 (density smaller than van Hove 
filling), deserves to be noticed. In the first case the two 
graphs on the left in Fig. show that the relative reduction 
of the local susceptibilities is stronger for Fermi surface 
points near to XM (9 = 9 m i n ), compared to Fermi sur- 
face points closer to the diagonal TM (9 w 7r/8, 6 = 7r/4). 
The situation changes if the calculations are done exclud- 
ing umklapp processes of all kinds from the second order 
contributions, see the two graphs on the right in Fig. [?|. 
This last fact already indicates the importance of elastic 
umklapp scattering at the Fermi surface that is allowed 
for densities larger than the van Hove density. The sec- 
ond case is when fi < —At 1 and elastic umklapp scattering 
does not enter the low energy physics. The two graphs on 
the left in Fig. || show a qualitative behavior of the lo- 
cal spin susceptibility that is similar to that observed for 
fi > —At' . But the strongest reduction of the local charge 
susceptibility is now shifted to 9 m tt/8 and 9 = tt/A and 
the role of the umklapp processes becomes marginal at 
these densities. 
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Fig. 8. Same as in Fig. |^ for band filling less than the van 
Hove filling such that elastic umklapp processes at the Fermi 
surface are absent. 



5 Ladder summations 

The results in the previous sections were obtained with 
the Landau interaction function calculated to second or- 
der in U. Although we observed a general suppression of 
the charge and spin susceptibility with respect to the non- 
interacting values, the effects were weak and no clear ten- 
dencies towards partial incompressibility could be found. 
In this section we extend the perturbative calculation of 
the Landau function to summations of infinite numbers 
of diagrams. Our motivation for doing this is twofold: a) 
A severe defect of the second order approach is that the 
contribution of the s-wave pairing channel to the interac- 
tion between the quasiparticles is typically overestimated 
to this order. In fact f a (kF,k' F ) = — U 2 xppO*-f + k^,) 



diverges at T — when 



-kp due to the singu- 



lar behavior of the Cooper graph in systems with parity 
e(k) = e(— k), that is clearly visible from the Fig. [|. In the 
following we will remedy this defect by summing up the 
ladder diagrams generated by the particle-particle graphs, 
the result is shown in the Fig. ^. In this way for U > 
the net contribution of the Cooper channel at T = will 
be zero, b) Furthermore we know that the Hubbard model 
close to half filling is close to an antifcrromagnetic insta- 
bility. This effect can be modeled in a perturbative way 
by summing up an infinite number of bubble and ladder 
diagrams in the particle-hole channel, as described below. 
By this we can enhance the proximity to the AF insta- 
bility and simultaneously study whether the local charge 
and spin susceptibilities exhibit significant changes. We 
should add that our calculation does not aim at a quan- 
titative determination of the phase diagram of the model. 
Therefore we do not worry too much about the question 
whether we neglect certain contributions by selecting this 
particular class of diagrams. Our goal is rather to include 



f(Q, 9') 








f\e, e') 




Fig. 9. Numerical solution of the symmetric f s (0, 0') and anti- 
symmetric part f a (8, 8') of the Landau interaction obtained by 
the ladder and bubble summation for t' = 0.3t, fi = — 3.83t', 
U — 1.5t and T = 0.04t. The unphysical divergences of 
f a ' s {kF,— kf) found with the second order approximation in 
Fig. are suppressed by the ladder summation. 



and enhance a specific type of scattering processes and to 
determine the effects of these processes on the Fermi liq- 
uid interaction and derived quantities. In the same spirit 
we neglect possible self energy corrections in these calcu- 
lations. Within the Fermi liquid framework this concerns 
mainly the effective mass of the quasiparticles, but a con- 
sistent treatment should also take into account the Fermi 
surface shift caused by the interactions and lifetime effects 
on the virtual excitations in the ladder diagrams. However, 
since generally in spatial dimensions higher than one the 
self energy does not develop singularities without clear 
signs in the interactions and since we seek for qualitative 
results that do not depend on details of the model, we 
ignore the selfenergy effects and focus on the interaction 
function. This has the advantage that the analysis remains 
relatively simple and lucid. 

The graphs contributing to the particle-particle and 
particle-hole ladder and bubble summations are shown in 
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Fig. 10. The only graph contributing to f^^(k.p, lt' F ) to 
the second order is the one-loop particle-hole bubble. This 
implies that the ladder summation for f^^ involves only 
odd powers of the bubble graph. 
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Fig. 10. Graphs contributing to the Landau interaction func- 
tions generated by the ladder summations of the first and sec- 
ond order graphs. 



Thus the Landau interaction functions can be expressed 



as 



/t,t 

/t,J. 



1-U*Xph^f-U f Y 

U U 2 xph^f-^ f ) 



(19) 



1 - UxppO^f + k' F ) 1 - UxphO^f - k' F 



(20) 



With the Landau interaction function obtained above Eqns. 
(|l|) and @ are solved again. Contrary to the second 
order approximation the local values of the spin suscepti- 
bility diverge to oo for U > U c with a critical U c (n) indi- 
cating a ferromagnetic instability of the system near the 
van Hove filling. This increased tendency towards ferro- 
magnetism can be explained by the suppression of singlet 
pairing by the ladder summation. In fact the same un- 
physical divergence of the second order particle-particle 
diagram was the main reason for the suppression of the 
spin susceptibility in the second order results. However, 
a more sophisticated ladder summation using a renormal- 
ized pair scattering jnj or an renormalization group ap- 
proach |],[l2j will tend to suppress the spin susceptibility 
again. This is because singlet pairing tendencies will be 
generated in higher angular momentum channels, espe- 
cially in the d,^^ -channel. 

The behavior of the local charge susceptibility when 
the ladder summation is included is a consistent exten- 
sion of the second order calculation with a stronger sup- 
pression of Xc(kp) near the saddle points. The graphs on 



the left in Fig. |TT] show Xc(kp) with close to a sad- 
dle point for four different values of U. At U = 1.62t ( 
thin dashed line), the local charge susceptibility goes to 
zero. This does not occur if the umklapp processes are ex- 
cluded from the perturbation expansion, as shown in the 
right panel of Fig. 
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The comparison between these two 
last graphs confirms the fundamental role of umklapp pro- 
cesses for the suppression of the charge susceptibility near 
the saddle points. In that sense the ladder summation re- 
sults for the Landau function lead to a similar behavior of 
the charge compressibility as the one observed in one-loop 
renormalization group treatments. 




X„ P 



0.05 




0.05 



Fig. 11. Left: k-space local values of XcQ^f) for k^ near a 
saddle point for different values of U: U/t = 0.5 thick line, 
U/t = 1.2 thick dashed line, U/t = 1.91 thin line, U/t = 2.05 
thin dashed line. Right: results obtained excluding umklapp 
processes from the perturbation expansion. For all data t' — 
0.3, fx = -3.99*' and T/t = 0.04. 



Another similarity with the renormalization group cal- 
culations are the tendencies towards deformations of the 
Fermi surface, also known as Pomeranchuk instabilities 
and originally suggested for the t-t' Hubbard model by 
Halboth and Metznerjl^] . Such an instability spontaneously 
breaks the fourfold symmetry of the electronic dispersion. 
Close to the van Hove filling the main energy gain comes 
from lowering the kinetic energy by pushing one saddle 
point deeper below the Fermi level. The Landau energy 
functional for this process can be written, as quadratic 
form in the Fermi surface shifts <5s(kp), 



SE = 



1 



(2tt)< 



dl F L s (5s)6s(k F ) 



(21) 



where L s is defined by ji~3|) . Negative eigenvalues of the 
linear integral operator L s imply a lowering of the energy 
( pT| ) by a suitable deformation of the Fermi surface. For 
the specific case t'/t = 0.3, fi/t' = -3.99 and T/t = 0.04 a 
negative eigenvalue appears already at U/t = 1.91, slightly 
before xA^-f) is suppressed down to zero close to the sad- 
dle points at U/t = 2.05. The corresponding eigenvector 
is shown in the graph to the left of Fig. 12 and signals a 



deformation 6s(9) of the Fermi surface breaking the point 
group symmetry of the square lattice, as shown schemati- 
cally in the plot on the right of the same figure. No attempt 
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was made to estimate the actual magnitude of the Fermi 
surface deformation. 




Fig. 12. Left: unstable eigenvector 8s (6) of the operator T s -i 
with eigenvalue A = — O.Olt < obtained for t' = 0.3t, /j, — 
— 3.99t', and U/t = 1.91. Right: Pomeranchuk deformation of 
the Fermi surface cx v(6)8s(9). 



In the case t' — 0.3t, the above scenario is observed 
for a very large range of chemical potential values, i.e. be- 
tween fi/t' = —5.66 and fi/t' = —3.33. Note however that 
this instability is only one example out of many potential 
instabilities of the 2D Hubbard model close to half filling. 
A direct comparison of several types of instabilities can 
be found in Ref. ||. 

6 Current carried by the quasiparticle and the 
Drude weight 

According to Landau-Fermi liquid theory the current Jk F 
carried by the quasiparticle kp is given by fill 

/dV 
^-r(k F ,k^)v k , F . (22) 

The first term of (0) describes the current of the bare 
quasiparticle and the second one the "backflow" of the 
medium around it. The Drude weight is connected to the 
quasiparticle current via |h| 

D =7Te J^; v ^ J ^ (23) 

With a simple symmetry argument it is easy to show that 
the off-diagonal components vanish. This remains true 
even when a Pomeranchuk instability breaks the point 
group symmetry of the square lattice, inducing a tetrago- 
nal symmetry. The Pomeranchuk instability merely mod- 
ifies the diagonal elements D aa and D bb , which would 
take different values. In this section we will assume that 

no Pomeranchuk instability affects the system, so that 
D aa =D bb = D 

Note that in a Galilean invariant system the total cur- 
rent is a constant of motion and does not change when 
the interaction is switched on adiabatically. But as a lat- 
tice model the Hubbard model does not belong to this 



class of systems, thus the magnitude of D is related to the 
strength of the interaction U. Fig. [l3| shows the change 
of the Drude weight when the Landau interaction func- 
tion that is used to solve (|2^) and ( p3[) is approximated 
by the ladder and bubble summations (|19|) and (pp|). The 
thick lines show D(U) / (e 2 Tr) for three different values of 
the chemical potential fi. With increasing interaction U, 
the Drude weight is reduced and goes to zero for a fi- 
nite U . The reduction occurs for smaller interactions the 
closer the electron density is to the van Hove filling. We 

0.5 r 

D/(e 2 ;i) 

0.4 -- - 



0.3 - 



0.2 




Fig. 13. The thick lines show D(U) obtained by the ladder 
summation for three different values of fi: [i/t' = —3.33 dashed 
line, fi/t' — —3.99 dash-dotted line, fi/t' = —5.66 continuous 
line and T=0.04t. The thin lines show the corresponding results 
obtained excluding the umklapp process. 



note that in our calculation the Drude weight vanishes 
only at interaction strengths that are larger than the crit- 
ical value necessary to induce the Pomeranchuk instabil- 
ity (U p : U/t = 2.075 for fj, = -3.33i, U/t = 1.91 for 
fj, = -3.99f and U/t = 2.6 for /i = -5.66t). Thus the 
above results have no clear physical meaning for U > U p . 
The behavior of the Drude weight for U < U p , however, 
provides some useful information: as shown in Fig. [l3], its 
decrease is strongest at the van Hove filling, [i = —At'. 
In addition, the results confirm again the important role 
played by the umklapp processes. In fact, if they are ex- 
cluded from the calculation (thin lines of the Fig. |l3|), the 
magnitude of the Drude weight even increases with the 
strength of U . 

It is interesting to determine the scattering processes 
responsible for the decrease of D. The continuous lines 
plotted in Fig. show the symmetric Landau interaction 
function f(9, 6*"for fj,/t' = —3.33 that has been used to 
solve (p2[) for different values of 9 and as function of 6'. 
The thick vertical lines mark the different values of 8 and 
the thin lines are plotted to emphasize the relevant peaks 
of f(9, 9') for the given 8. The insets shows the locations of 
kp(0) (thick line) and k' F (9') (thin lines) where the peaks 
in f(6, 9') occur. 

Fig.|l4| shows that the peaks in the Landau function 
for fillings larger than the van Hove filling are tied to the 
Fermi surface points ki?, k^ which can be connected by 
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f s (e,e') 




Fig. 14. The continuous lines show f(9,9') for fi/t' = —3.33, 
U — 2.075t and T = 0.04t, obtained for three different values 
of 9, as a function of 9'. The thick vertical lines mark the 
different values of 9 and the thin vertical ones are plotted to 
emphasize the relevant peaks of f(9,9'). The insets show the 
locations of the Fermi surface points marked by the thick and 
thin lines. The dashed lines show the corresponding results of 
f(9,9') obtained excluding the umklapp process. 



wavevectors close to (tt,tt), i.e. qo = kp — k'p ~ (±7r, ±7r). 
These processes are enhanced in the ladder and bubble 
summation of the terms involving xph(1o) (@)- For fillings 
smaller than the van Hove filling and —6.6 < fi/t' < —4 
the Fermi surface is partially nested and the peaks appear 
when q = kp — k' p is equal to the corresponding nesting 
vectors (that are smaller than (it, 7t)). For fj,/t' < —6.6 
the Fermi surface is strictly convex and the Drude weight 
increases with the strength of U. 

Without umklapp processes the relevant peaks of f s 
disappear, the corresponding Landau interaction (thin dash- 
ed lines of the Fig. 14) causes a "backflow" of the medium 
in the direction of the bare quasiparticle current implying 
the increase of Jk F for anyvalue of kp. Thus the Drude 
weight will increase. 



7 Conclusions 

We have presented a detailed study of the Landau interac- 
tion function /(k, k') and wave-vector resolved compress- 
ibilities within Fermi liquid theory for the two-dimensional 
Hubbard model. This model is often considered as mini- 
mal model for the description of high-T c superconductivity 
in the copper-oxide materials. 

A less surprising conclusion of this work is that a Fermi 
liquid picture alone can hardly account for the observed 
anomalies that occur in the underdoped high-T c materials 
when the doping is reduced towards half-filling: all effects 
are rather weak and one has to resort to stronger inter- 
actions U or move closer to a magnetic instability (e.g. 



by ladder summations) to obtain sizable effects such as 
a significant reduction of the compressibility on parts of 
the Fermi surface. No rapid decrease of the spin suscepti- 
bility reminiscent of spin gap formationjll) can be found. 
Similar conclusions regarding the high-T c cuprates were 
reached by Kim and Coffey|l™ who analyzed the spectral 
function of the 2D Hubbard model within random phase 
approximation. 

On the other hand, our weak coupling analysis pro- 
vides insights into the physical processes which lead to 
deviations from the conventional picture. We find that 
many of the observed tendencies towards partially insu- 
lating behavior - as we believe an essential ingredient of 
the underdoped cuprate superconductors - are related to 
elastic umklapp processes. These processes start to affect 
the low energy physics upon reducing the doping when the 
Fermi surface touches the Brillouin zone boundaries. As 
we know from earlier renormalization group studies |^,^| of 
the same model, strong umklapp processes with momen- 
tum transfer (-7T, it) between the saddle point regions lead 
to an interesting mutual reinforcement of antiferromag- 
nctic and rf-wave pairing correlations. Furthermore the 
renormalization group treatment yielded a regime with 
a strong suppression of the charge compressibility near 
the saddle points. A similar effect was observed by Ot- 
suka et al.jlTj] in quantum Monte Carlo calculations for 
U = At. In our Fermi liquid analysis we clearly identify 
the elastic umklapp processes as being responsible for a 
strong reduction of the Drude weight and the angle re- 
solved charge compressibility on parts of the Fermi sur- 
face. In one-dimensional systems it is well known that 
elastic umklapp processes imply these effects. Our calcu- 
lations shows that in the 2D Hubbard model there are 
observable tendencies in the same direction. Although our 
perturbative analysis does not allow to conclude that umk- 
lapp processes are the ultimate cause for the pseudogap 
features at stronger interactions, they will at least con- 
tribute to a certain extent. Furthermore there is no obvi- 
ous reason why the importance of the umklapp scattering 
should decrease again for stronger interactions. 

Our perturbative calculations show tendencies which 
can be interpreted as precursors to stronger reductions 
in the compressibility and Drude weight at stronger in- 
teractions, features which appear in some phenomenolog- 
ical descriptions of the socalled pseudogap region of the 
phase diagram of the underdoped cuprates. In this way 
they support the proposal that the critical doping of the 
cuprates that separates the underdoped and overdoped 
regimes 0, is determined by the appearance of elastic 
umklapp scattering with momentum transfer (tt, it) con- 
necting the Fermi surface regions near to the saddle points. 
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port by the German Science Foundation (DFG). Compu- 
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